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Abstract. In this paper we study a generalization of the Fibonacci sequence 
in which rabbits are mortal and take more that two months to become mature. 
In particular we give a general recurrence relation for these sequences (improv- 
ing the work in [6]) and we calculate explicitly their general term (extending 
the work in [7])- In passing, and as a technical requirement, we also study 
the behavior of the positive real roots of the characteristic polynomial of the 
considered sequences. 
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1. Introduction 

Fibonacci numbers arose in the answer to a problem proposed by Leonardo de 
Pisa who asked for the number of rabbits at the n*'' month if there is one pair of 
rabbits at the O*'' month which becomes mature one month later and that breeds 
another pair in each of the succeeding months, and if these new pairs breed again 
in the second month following birth. It can be easily proved by induction that the 
number of pairs of rabbits at the n*'' month is given by /„, with /„ satisfying the 
recurrence relation: 

/o = /i = 1. 

fn = fn-1 + fn-2, Vn > 2. 

It is not the point here to state any of the many properties of these numbers (see 
[8] for a good account of them), nevertheless we will recall that if ri < r2 are the 



roots of the polynomial g{x) = a 

(1) fn 



1 then we can see that: 



n ~r2 r2 - ri 



In [7] the fc-generalized Fibonacci numbers /„ are defined as follows: 

fi^^ = 1, VO < n < fc - 1. 



fe 



i=l 



In this paper Miles proves, among other results, that if ri, . . . , are the (distinct) 
roots of gk{x) = x'^ — x'^~^ — ■ ■ ■ — x — 1 then: 



(2) 



\ 



1=1 



n (n-r,r' 



\l<i<fc 



which, of course reduces to ([T|) if we set k = 2. 
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Later, in [61, Hoggat and Lind consider the so called "dying rabbit problem", 
previously introduced in [1] and studied in [2] or , which consists in letting rabbits 
die. In the mentioned paper Hoggat and Lind consider a pair of rabbits at the 0*'* 
time point which produces _B„ pairs at the n*'* time point and which die at the m*'' 
time point after birth. Then if T„ is the total number of live pairs of rabbits at the 
n*'* time point, and defining B{x) = ^ Bnx'\ D{x) = x"' and T{x) = ^ r„a;", 

n>0 n>0 

they prove that 

^ > il-x){l~B{x))- ^ 
and use this formula to find recurrence relations for T„ in various cases. 

The goal of this paper is to give a new look at the dying rabbit problem. In the 
second section we study a family of polynomials, focusing on the behavior of their 
positive roots. Although motivated by technical requirements, this study turns out 
to be of intrinsic interest. In the third section we will find a general recurrence 
relation for the sequence arising in this problem (which is only given in ^ for some 
particular cases) and we will deduce an explicit formula (which also generalizes the 
work by Miles) for the total number of live pairs at the n*'' time point. Finally, in 
an appendix, we give a procedure written using Maple® to calculate terms of the 
considered sequences. 

2. A FAMILY OF POLYNOMIALS AND THEIR ROOTS 

Given natural numbers /i, fc > 1 we define the following polynomial: 

gkAx) = x''+''-^ - x"-^ ^-1. 

In this section we will study, in some sense, the behavior of the roots of this polyno- 
mial in terms of k and h. In particular we will be interested in the unique positive 
real root of gkji{x). 

Proposition. gk.h{x) has a unique positive real root ak,h which lies in (1,2). 

Proof. Just apply Descartes' rule of signs together with Bolzano's theorem. □ 

Remark. Note that ai./j = 1 for all h > 1. 

As a consequence of the previous proposition we have the following lemma. The 
proof is elementary and we omit. 

Lemma. Let h,k>l and < ?; G M, then y > ak.h if and only if gk,hiy) > 0. 

Now, for every fixed h we consider the sequence {ak,h}k>i- Similarly, for ev- 
ery fixed k we have a sequence {afc^;i};i>i. The following proposition studies the 
monotony of these sequences. 

Proposition. Let h,k>l. Then: 

(1) ak,h < ak+i,h- 

(2) ak,h > CKk^h+i, md the equality holds if and only if k = 1. 

Proof. (1) By definition we know that gk+i,h{(Xk+i,h) = 0. Now, we have that 

k+h-i "fc+i\ "fc+i.fe + Q^fc;Lfe + --- + «fc+i,/» + i ^ k-i ^ ^ 

o^k+i.h ^ > "fe+i,h + ■ • • + 

ak+i,h afc+i,/i 

Q.k+i,h + 1, so gkJi{ak+iJi) > and the result follows from the previous 
lemma. 
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(2) Again by definition gk,h{oik,h) — and we can write gk,h+i{ock,h) 



a 



k-l 

k,h 



k+h 
'^k,h 



■ - ak,H - 1 = 4%^ - 4+!^-' = 4%''-\^k,h - 1) > 0, with 
the cquahty holding if and only if = 1. An application of the lemma 
completes the proof. 

□ 

Now, as we have seen that both sequences {oik,h}k>i and {(Xk,h}h>i are mono- 
tonic, and as they arc clearly bounded, they must be convergent. The next result 

is devoted to calculate their limits. 

Proposition. (1) lim ak,h — cth for all h > 1, where ah is the unique positive 

k — ^oo 

root of the polynomial Ph{x) = — x^~^ — 1. 
(2) lim ak^h = 1 for all k > 1. 



Proof. (1) Let us fix /i > 1. Then for any A; > 2 we have 



k+h-l 



1 + o:k,h + 



k-l 
^k,h 



a 



k.h 



- 1 



ak,h - 1 



and thus a 



k.h 



^k,h 



1 



Now, as wc know 



that ah = lim ak^h > 1 it is enough to take limits in the previous equality 

k — *oo ' 

to obtain the result. 
(2) Let us fix now fc > 1. Then for any h > 2 we have a^^''"^ 

.k-l 



1 + ak,h + 



+ a^ so, we take logarithms in both sides to obtain the equality 



log(l + afe,?, + --- + 4_^i) 



logQ;fc,/t = 



Finally, if we call /Jfe = lim ak^h 



k + h-l 

and take limits in the previous expression we arrive at log = for every 
fc > 1 and the proof is complete. 

□ 

The previous propositions can be summarized in the following diagram: 



Q!l,l 


< 


Q!2,l 


< 


03,1 


< 


Q!4,l 


< . . 


. < 


Ofc,! 


< . . 




Oil 


11 




V 




V 




V 






V 






V 


ai,2 


< 


0:2,2 


< 


03,2 


< 


0:4,2 


< . . 


. < 


Ofe,2 


< . . 




02 


II 




V 




V 




V 






V 






V 


"1,3 


< 


^2,3 


< 


03,3 


< 


0:4,3 


< . . 


. < 


o:fe,3 


< . . 




0:3 


II 




V 




V 




V 






V 






V 


Q!l,4 


< 


Q!2,4 


< 


03,4 


< 


04,4 


< . . 


. < 


afc,4 


< . . 




0:4 


II 




V 




V 




V 






V 






V 




< 


0:2,5 


< 


03,5 


< 


04,5 


< . . 


. < 




< . . 




05 


II 




V 




V 




V 






V 






V 


1 




V 




V 




V 






V 






V 


ai,h 


< 


o:2,h 


< 


o:3,h 


< 


04,/» 


< . . 


. < 


Olk,h 


< . . 
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II 




V 




V 
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V 
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1 
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1 






1 



Where ah is the unique positive root oiph{x) 
is strict. 



x'^-x'^-^- 



1 and every inequality 
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Before we go on, we introduce a result by Cauchy (see [3]) which wiU be useful 
in the forthcoming. This result gives a bound on the modulus of the roots of a 
polynomial with complex coefficients. We present it without proof. 

Theorem. Let f{z) = z" + a„-iz"^^ + ■ ■ • + aiz + ao be a polynomial with complex 
coefficients (oi ^ for at least one i) and let Z[f(z)] the set of its complex roots. 
Then Z[f{z)] C {z E C \ l^l < 77} where rj is the unique positive root of the real 
polynomial f{x) — x" ~ \an-i\x"^^ — • • • — |ai|a; — |ao|. 

Corollary. Z[gk,h{x)] C {z G C | \z\ < ak,h}- 

Proof. Since gk,h{x) = gk.hix), it is enough to apply the previous theorem. □ 

Now, we can refine the previous corollary in the following way. 

Proposition. If gk^h(w) = then \w\ < ak,h o,nd moreover, \w\ ~ oikji if and only 
if w = ak.h- In other words, ak,h is the largest root of gk.h{x). 

Proof. In what follows we will put a — ak,h. Let w be a root of gk.h{x) such that 
— a. As gk,h{w) = we have that w^~^{'w'^ — I) = + - • ■+w+l and, in par- 
ticular, jw*^""'^ pju''' — — Ip = — Ip. If we put w — a cosO+ia sin 9 we ob- 
tain the equality a^^^"^ {{a^ cos hO - \ f + a^^' sin^ hO) ((a cos - if + o? sin^ B) = 
{(a^^^ cos{k ~ 1)6 ~ if + o?^^'^ sin^(A; - 1)6') and operating we get a^^^'^(a^^ -I- 
1 - coshd){a^ + 1 - 2a cos 6*) = a^^-"^ + 1 - 2a'=-i cos(fc - 1)6*. 
Let us now define polynomials p{x) and q(x) by: 

then, it is easy to see that coshd — p{cosd) and cos(A; — 1)6' — q{cosd). Moreover, 
p{l) — q{l) = 1 and if we define r{x) to be 

r{x) = a2fe-2(c,2'i + 1 _ 2a''p{x)){a^ + 1 - 2ax) - (a^^-a + 1 _ 2a''-^q{x)) 

then, gk^h{a) — implies r(l) = 0. Also we can see that r{x) has no roots in 
(—1, 1). Thus, if ui e C is a root of gk.h{x) with \w\ = a it must be real and as 
gkjii—oi) ^ 0, we conclude that 

Z[gk,h{^)] \ {akji} C {z e C I |z| < ak,h} 
and the proof is complete. □ 



We will finish this section with the following proposition which will be of great 
technical importance in the next section. 

Proposition. All the roots of gk.h o,i"c distinct. 

Proof. We will show that gk,h{x) and g'^. j^{x) have no common root. To see so it is 
enough to prove that g.c.A.{gk,h,g'k h) ^ constant. But we can use Euclid's algo- 
rithm repeatedly (multiplying by appropriate constants in every step if necessary) 
to arrive to the result. □ 



THE DYING RABBIT PROBLEM REVISITED 



5 



3. The dying rabbit sequence 



As we mentioned in the introduction, we are interested in generalizing the Fi- 
bonacci sequence by considering that rabbits become mature h months after their 
birth and that they die k months after their matureness. We will denote by C, 
the number of couples of rabbits at the n*'* month. Obviously we have: 



ik,h) 



1. 



Now let us denote by cl!^^ the recurrence sequence defined by: 



c. 



(h) 



ci%Wn>h 



then an easy induction argument let us see that; 



r'ik,h) _ 



C\ 



{h) 



^{k,h) ^(k,h) ^ik,h) 
^^n-h "1 '^n-k-h+1' 



h- 
2 



if < n < A: 
if n> k + h- 

Example. (See the appendix) 

• If A; = 3 and h = 2, the beginning terms of Cn'^^ are: 

1,1,2,3,4,6,9,13, 19,28,41,... 

• If k = 7 and h = 4, then the beginning terms of Cn'^^ are: 

1, 1, 1, 1, 2, 3, 4, 5, 7, 10, 13, 17, 23, 32, . . . 

Remark. The characteristic polynomial of the recurrence sequence C^'^^ is pre- 
cisely the polynomial gk,h{x) studied in the previous section. 

If we denote by ri, r2, . . . , rfc+h_i the (distinct) complex roots of gk,h{x) it follows 
from section 2 and from well-known facts from the theory of recurrence sequences 
that there exist constants ai, 02, ... , ak+h-i such that: 

^(^■^^ = a^r^ + a^r^ + ■■■ + au+h-irl+h-i 
where we can suppose that ri = a^^h- In particular we can calculate those constants 
solving the system of linear equations given by: 

k+h-l 

(3) ^ air\ = c\^\ 0<l<k + h-2. 

which can be expressed matricially: 
/I 1 ... 



r-2 



k+h-2 „k+h-2 



1 \ / 

Tk+h-l 

Ji+h-2 



ai 
a2 



(h) 



(h) 

k+h-2 



and which has unique solution because all the arc distinct. 

To solve this system of equations we will use Cramer's rule. Recall that if we 
put 

1 1 •■• 1 

ri r2 ■■■ Vk+h-i 



V 



^k+h-2 „k+h-2 
1 ^2 



r.k+h-2 
'^k+h-1 



n i'-^j) 

k+h-l>i>j>l 
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Dn = 



1 

„k+h~2 



1 

Tn-l 
'n-1 



Co 



1 

rn+1 

k+h-l 



1 

rk+h-1 



k+h-2 ''„+! 



r 



k+h-2 
k+h-1 



then a„ = . So it is enough to find Z3„ for n — 1, . . . , k + h — 1. We will work 
out the case n — 1 completely, the other cases being analogous. Also note that 
we have replaced the values Cj''"' {j — 0, . . . , k + h — 2) hy arbitrary constants Cj 
{j — 0, . . . , k + h — 2) , a,s their value is of no importance when solving ([3]) formally. 

To calculate Di we first need the following generalization of Vandermonde de- 
terminant which can be found in [5J (Lemma 2.1). 

Lemma. //e„ is the n*"^ elementary symmetric polynomial in the variables {xi, . . . ,x 
then: 

1 ... 1 1 



J]^ (xi-Xj) e„_i(xi,...,a;„). 



. n>i>j>l 



Now, if we apply this lemma and we expand the determinant Di by its first 
column we obtain: 

k+h-2 I \ 

-Di = ^ {-^)^Ci\ (ri -Tj) efc+,,_2-i('^2,-.-,rfe+/i-i) 

;=0 \k^h-\>i>j>2 j 

and, consequently: 

^ ^ k+h-2 

(4) ai = -r^ = } : ^ (-l)'Ciefc+,,_2-i(r-2, . . . ,rfe+ft_i). 

I _[ (ri — ri) 

k+h-l>i>2 

So, we are now interested in calculating the values ej(r2, ■ ■ ■ ,rk+h-i) for < 
j < k + h — 2. From Cardano's formulae and considering that ri, . . . , rj;^/i_i are 
the roots of gk,h{x) we have that: 

eo(ri, . . .,rk+h-i) 1- 
ei(ri, . . . ,rk+h-i) = ■■■ = e/,._i(ri, . . . ,rk+h^i) = 0. 
e,(ri, . . . , rk+h-i) = (-1)'+' V/i < s < fc + - L 
On the other hand, the following lemma is easy to prove. 

Lemma. et(x2, . . . , x„) = ^(-l^+i^^i^^^i^^^^^ VO < t < n. 

We put this together to obtain: 

60(^2, • • • = 1- 

k ^ 

e.(r2, . . . , rfc+ft_i) = (-1)^ ^ VI < s < - 1. 
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k-\-h—l — s 

es{r2,...,rk+h-i) = {-iy T ^h<s<k + h~2. 



and summing the geometric series: 

eo(?'2, • ■ - ^rk+h-i) = 1 

r'[ - 1 

rJ+''-i-^(ri - 1) 



e,(r2, . . . , n+h-i) = (-1)' k+h-i-s, VI < s < - 1. 



^k+h-l-s _ ^ 

es{r2, ■ ■ . = (-1)^ k+h-i-s/ Vh<s<k + h-2. 



Finally if we substitute in ^ we get 



ai 



— 1 „ ri — 1 



k+h-l>i>2 



Reasoning in a similar way and taking into account the symmetry of the we 
can calculate a„ for every l<n<k + h— \. In fact: 



(-If 



k+h+n-l 



T-r T-r ~ y^^'^TTTT 7+ Q " + C'fc+/i_2 



Remark. Observe that the previous considerations are only valid for k,h > 2. For 
the case /i = 1 we refer to Miles' paper [7] and the case fc = 1 is trivial. 



that \ri\ < \ri\ for all z > 2 we have that lim "^^.^ — fi — ak^h- This generalizes 



Remark. It is interesting to observe that ai 7^ 0. As a consequence and recalling 

n— >oo r^i 
f 1 

the fact that = $ where /„ is the n*^ Fibonacci number and $ is de golden 

Jn 

section (note that 02,1 = $). The previous expression for a„ also generalizes the 
work by Miles just by setting h = 1. 

Example. (Padovan sequence) 

Recall that the so-called Padovan sequence is defined by 

Po = Pi=P2 = 1. 

Pn = Pn-2 + Pn-3, Vn > 3. 

(2 2). 

Thus, it is clear that in our notation P„ = Cn ' with the initial conditions Cq = 
Ci ^ C2 — I- So we can apply our previous results to obtain that: 

^ ^ rl+r, + l r|+r2 + l^„ + + I 

2ri + 3 1 2r2 + 3 ^ 2r3 + 3 ^' 
which was already known to hold. 

If we keep the same recurrence relation but replace the initial conditions by Pq = 
3, Pi = 0, P2 = 2 we obtain the so-called Perrin sequence, whose general term can 
be again calculated with our formulas to obtain: 
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Finally, if we keep our original initial conditions, that is, Cg =(7^ =1,(72 = 
2, then the general term of our Padovan-Perrin like sequence turns out to be: 

^" ~ 2ri + 3 1 + 2r2 + 3 ^ + 2r3 + 3 ' 
4. Appendix 

In this appendix we give a short and easy procedure, written with Maple®, which 
calculates any number of terms of Cn''^^- It goes as follows: 
dr :=proc(k,h,t) 
local i; 

for i from by 1 to h-1 do c(i):=l end do; 

for i from h by 1 to k+h-2 do c(i) :=c(i-l)+c(i-h) ; end do; 

for i from k+h-1 to t do c(i) :=siim(c(n) , n=i-k-h+l . . i-h) ; end do; 

print (seq(c (n) ,n=0. .t)) ; 

end proc: 
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